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The distribution of temperature and emission measure in the stationary heated 
solar atmosphere was obtained for the limiting cases of slow and fast heating, when 
either the gas pressure or the concentration are constant throughout the layer depth. 
Under these conditions the temperature distribution with depth is determined by 
radiation loss and thermal conductivity. It is shown that both in the case of slow 
heating and of impulsive heating, temperatures are distributed in such a way that 
classical collisional heat conduction is valid in the chromosphere-corona transition 
region of the solar atmosphere. 
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I. INTRODUCTION 

Recently appeared a number of papers [1],[2] claiming that the chromosphere- corona 
transition region of the solar atmosphere should be considered in the non-collisional approx- 
imation. In these papers it is said that ion-acoustic turbulence is the reason of differentiation 
of solar plasma into two regions with high (T e £ 10 6 K) and low (T e ^ 10 4 K) electron tem- 
perature, correspondingly. 

In this paper we present the solution of the equation of balance between the thermal 
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heating and radiative cooling with classical electron conductivity. This solution explain dif- 
ferentiation of solar plasma to the high and low electron temperature. The characteristic 
thickness of the chromosphere-corona transition region is greater then thickness correspond- 
ing to the free path for thermal electron collisions and such temperature distribution agrees 
with observed solar radiation. 



II. CALCULATION OF TEMPERATURE DISTRIBUTION WITH DEPTH 

In this paragraph we obtain temperature T distribution with depth £ for the heating 
solar plasma by thermal flux, taking radiation loss into account. 
The thickness is defined as 



n(x)dx, (1) 

where n(x) is the number density of particles. The thickness represents the number of atoms 
in the column of unit cross-section along the direction of thermal flux. The distribution 
T = T(£) can be obtained with the help of the equation of balance between the thermal 
heating and radiative cooling: 

d ( dT\ T /rrtS 2 



n T ( { Kn dz) = L{T)n P ~' (2) 

where k is electron thermal conductivity, 

. /2 1.84 x i(T 5 T 5 / 2 
^oT 5 / 2 = — , (3) 

[3], [4], L = L(T) is the distribution of radiative energy losses with temperature (see Fig.l), 
P = LiT^n 2 - stationary thermal heating at the infinity. 

Parameter L of the equation (2) are known: L = L(T) was taken from [5]. To determine 
the number density dependence on the temperature n = n(T) let us consider the cases of 
slow and fast heating. 

One has to specify doundary conditions for the temperature in the equation (2): 



dT 
~d£ 



= 0, T 
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According to [6], to obtain the dependence T = T(£) from the equation of the thermal 
balance (2) one has to multiply both parts of the equation (2) by k^. Then after simple 
manipulations we have the following expression for the thermal flux F: 



F = 2{L{T') - L(T oc ))K,n 2 dT'^j ' , 



(5) 



where the thermal flux F is defined as: 



^ dT 

F = - K „-. (6) 



To express (5) in dimensionless form we multiply (5) by 

e 2 

So 



-.1 



^OO K OO Too L ( Tqq ) 

where £oo, ^oo, T^, Koo are the units of the depth, concentration, temperature and electron 
conductivity, correspondingly. 
We obtain 

/ pT T 5/2 \ V2 

F=^ 2K 1 (L(T)-l)(—)n 2 dT) , (7) 
where In A is the Coulomb logarithm, 



In A 



In (1.24 x 10 4 (T Too) 3 / 2 /(n n^) 1 / 2 ) , T e < 5.8 x 10 5 K, 
In (9.44 x 10 6 (T T^/in n^) 1 ' 2 ) , T e > 5.8 x 10 5 K 



And K x 



K X = & L( ?<*>) . (8) 

Let us choose the units to be the following (recall that these values are simply our choosen 
units and not the actual values at the infinity): 

T^ = 10 A K, 

noo = 10 10 cm -3 , 

1 /2 

U = noo/oo = 3.15 x 10 15 cm" 2 , ^ = f r ^°° 2 ) = 3.15 x 10 5 cm. (9) 

\M i °oJ Tl oo/ 

From equation (6) we obtain the following expression for £: 

f Tc dT 

( = J T Kn Wr (10) 
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We choose the origin of ^-coordinate to be the point where the temperature has a fixed 
value T c , T c = 10 6 K ( f = when T = 10 6 K). 

Finally, to get the distribution T(£) we can find F(T) from equation (5) and then, after 
putting F(T) into (6), obtain £(T) and thus T(£). 

In the case of fast heating, when concentration is constant throughout the layer depth 
(n = const), let us set n = n^, or in dimensionless form n — 1. 

In the case of slow heating gas pressure is constant throughout the layer depth (p = const). 
In this case, because 

p = nk B T, (11) 



concentration becomes 



B = iSf (12) 



or in dimensionless form: 

n=\- (13) 



The dependence of the thermal flux on the temperature F = F(T) for the cases n = const 
and p = const is shown in Fig. 2. Here 

Foe = 7 = 425 erg/s. 

Soo 

The dependence of the temperature along the depth T = T(£) for the cases n = const and 
p = const is shown in Fig. 3. 

III. THE WIDTH OF THE TRANSITION REGION (<5£(T) AND ^ e 

COMPARISON) 

Classical collisional heat conduction [3], [4] is valid if two following conditions are satisfied: 

K<h= wk (14) 

where A e is the mean free path for thermal electron collisions, It is the characteristic scale 
of length on the temperature profile. 
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Condition (14) may be written as: 

^ > (15) 

where is the characteristic depth for equilibrium temperature distribution obtained above, 
and £ e is the thickness corresponding to the free path for thermal electron collisions. 

*< r > = &}• (16) 

k 2 T 2 

£e = n e Xe = -f— r- (17) 
7re 4 In A 

The dependence 5£ = 5£ (T) for the cases n = const and p = const and the dependence 
£ e = £ e (T) are shown in Fig. 4. 

As we can see from Fig. 4 5£ » £ e , that is the characteristic thickness at which tempera- 
ture is changing is more then thickness corresponding to the free path for thermal electron 
collisions in 400-500 times. For example, on the temperature T = 10 5 K temperature is 
changing on the 35 km and the free path for thermal electron collisions 70 m. 

Thus the chromosphere-corona transition region of the solar atmosphere shall be consid- 
ered in the collisional approximation. 

IV. STABILITY OF THE SOLUTION 

Linear theory of the thermal instability was constructed in monography [7] due to Field. 
The uniform medium in the thermal and mechanical balance in linear theory is characterised 
by 3 dimensionless parameters a, (3, 7. 

If the heating is such that energy gain per second to the gramm of substance is now de- 
pendent on temperature and density, and cooling of the medium becomes formed by volume 
radiation energy losses (n 2 L(T)), than the parameter a will depend only on temperature. 
And in this case a is the logarithmic derivative of L(T) 

, m s d InL 

a(T) = ItaT" (18) 
Parameter j3 characterises comparative significance of the thermal conductivity. If the con- 
ductivity is defined only by free electrons (3), then 
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/ 3( T) = l2_Jl!^VTL(T) (19) 

7 k B 

That is, /3 also depends only on the temperature; here \x is the effective molecular weight 
(for plasma with cosmic abundance of elements /i ~ l,44m#). 

Dependence of a and (5 on the temperature is shown in the Fig. 5. We assume 7 = 5/3 
in the equation (19). 

On figure 7 we can see regions where perturbations of the following types can be unstable: 
1) Isobaric perturbations, for which p = const. Regions where a < 1 correspond to isobaric 
perturbations. This mode is called condensation mode of thermal instability. 2) In regions 
where a < —3/2 adiabatic (entropy = const) perturbations are unstable. This mode of 
thermal instability is named wave or sonic. 3) In regions with a < isochoric (n = const) 
perturbations are unstable [8]. 

Let [7], 16: 

„ „l/2 1.3/2 T 3 /2 

e - — - 2 (20) 

Zt = y t = a ~ l Zp> ( 21 ) 

^ = -r = H P . (22) 

On the scales which are smaller than critical [7] ,26, thermal instability is stabilized by 
condactivity: 

£ cc = £p/3 1/2 (l-«r 1/2 , (23) 

U = U- 1/2 (-«-^rr 1/2 (24) 

7 - 1 

for the condensation and wave mode. These values of thickness which also depend only on 
temperature, are shown in the Fig.6. 

The values of characteristic thickness which correspond to the biggest increase rate of 
thermal instability are [7], 46: 

e„ c =(^ + ^)" /4 (u„)" 2 , (25) 



1^ i|-l/2/t t \l/2 
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(£, U) Vi (26) 
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for the condensation and wave mode, correspondingly; they are also shown in the Fig.6. 

The characteristic thicknesses corresponding to the stationary chromosphere heating in 
the cases p = const and n = const for the temperature profiles are shown in Fig.6. At all 
points of distribution the balance between the thermal heating and radiative cooling has its 
place. 

Let us look at Fig.6. If A < A cr , then perturbations are smoothened out by electron 
conductivity, if A ~ A cr , then pertrubations will grow. Note, that our solutions (in the cases 
= const and n = const) are crossed by the curve £ cc corresponding to condensation insta- 
bility. Then, if in the places of instability there is a spontaneous change of the temperature 
profile, the profile will return to its original position. 

So, we have undrestood that the temperature profile cannot be flatter than equilibrium 
plofile. It also cannot be steeper, because then the same temperature will be accumulated 
at smaller thickness, so the emission in this temperature range will be lower. 



V. THE TEMPERATURE PROFILE EMISSION 

The ability to emit of a certain region is named the emission measure (ME). 



ME = j n{ dl, (27) 
o 

here x is the length of emitting region along the line of sight, dl is the interval of length. 

Differential emission measure (DME) is the derivative of emission measure with respect 
to temperature. 

DME = i™. (28) 

Let us rewrite DME in terms of £: 

d ME nl dx n e d x d £ , . 

DME = ——— = = n e = n e -±- 29 

d T dT dT dT y ' 

Now, we can calculate the distribution DME = DME(T) for the cases n = const and 
p = const, using the dependences £(T) and n e = n(T). The result and measured DME 
points (for different lines) is shown in the Fig. 7. 
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VI. CONCLUSION 

The distribution of temperature with depth was found, assuming that the electon conduc- 
tivity had place, and at all points of distribution the balance between the thermal heating 
and radiative cooling had place. Our solution is stable (see IV) and observed UV-radiation 
can be explained by it (see V). 

The obtained results can be used to show that temperatures are distributed in such a 
way that classical collisional heat conduction is valid in the chromosphere- corona transition 
region of the solar atmosphere, because the characteristic thickness, at which temperature is 
changing greater then thickness, corresponds to the free path for thermal electron collisions. 
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FIG. 1: Dependence of the total radiative power loss on the tempetature L = L(T). 
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T, 10 K 

FIG. 2: Dependence of the thermal flux on the temperature F = F(T). 
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FIG. 3: Distribution of the temperature along the depth T = T(£). For the cases fast (n = const) 
and slow (p = const) heating. 



12 




FIG. 4: Characteristic values of the thickness versus the temperature for the equilibrium temper- 
ature distriution <5£ = <5£(T) for the cases fast (n = const) and slow (p = const) heating. And the 
thickness corresponding to the mean free path of thermal electrons £ e = (, e (T). 
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FIG. 6: The scales of condensation £ cc and wave £, cw perturbations, which is stabilizate by conduc- 
tivity ((23), (24)), and scales of perturbations with the highest rate of growth £, mc ,£, mw ((25), (26)) 
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FIG. 7: Dependence of the differential emission measure on the temperature DME = DME(T) 
and measured DME points (for different lines). 



